home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / dzasum.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.1 KB  |  49 lines
  1.       DOUBLE PRECISION FUNCTION DZASUM( N, ZX, INCX )
  2. *
  3. *     takes the sum of the absolute values.
  4. *     jack dongarra, 3/11/78.
  5. *
  6. *     .. Scalar Arguments ..
  7.       INTEGER                           INCX, N
  8. *     ..
  9. *     .. Array Arguments ..
  10.       COMPLEX*16                        ZX( 1 )
  11. *     ..
  12. *     .. Local Scalars ..
  13.       INTEGER                           I, IX, NINCX
  14.       DOUBLE PRECISION                  STEMP
  15. *     ..
  16. *     .. External Functions ..
  17.       DOUBLE PRECISION                  DCABS1
  18.       EXTERNAL                          DCABS1
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       DZASUM = 0.0D0
  23.       STEMP = 0.0D0
  24.       IF( N.LE.0 )
  25.      $   RETURN
  26.       IF( INCX.EQ.1 )
  27.      $   GO TO 20
  28. *
  29. *        code for increment not equal to 1
  30. *
  31.       IX = 1
  32.       IF( INCX.LT.0 )
  33.      $   IX = 1 - ( N-1 )*INCX
  34.       NINCX = IX + ( N-1 )*INCX
  35.       DO 10 I = IX, NINCX, INCX
  36.          STEMP = STEMP + DCABS1( ZX( IX ) )
  37.    10 CONTINUE
  38.       DZASUM = STEMP
  39.       RETURN
  40. *
  41. *        code for increment equal to 1
  42. *
  43.    20 DO 30 I = 1, N
  44.          STEMP = STEMP + DCABS1( ZX( I ) )
  45.    30 CONTINUE
  46.       DZASUM = STEMP
  47.       RETURN
  48.       END
  49.